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The solution of the time-dependent Schrodinger equation is discussed for a particle confined in 
half-space a: > with a linear potential V{x) — Kx in the following situations: (a) sudden removal 
of the wall and switching on the linear potential V{x) = Kx at 2: < 0, (b) sudden removal of the 
wall and the potential and (c) sudden removal of the potential. A brief discussion of the quantum 
statistic is presented. 
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I. INTRODUCTION 

The quantum mechanical efFccts resulting from time-dependent boundary conditions have been studied by some 
authors [J-ISi I6l4l3| . Diffraction in time was initially introduced by Moshinsky [l| by considering a situation involving 
a beam of particles impinging from the left on a totally absorbing shutter located at the origin which is suddenly 
turned off at an instant. In such an example, the transient probability current has a close mathematical resemblance 
with the intensity of light resulting from the Fresnel diffraction by a straight edge. See Q for a recent review. An 
interesting feature of the solutions for cutoff initial waves, occurring both in the free case [H and in the presence 
of a potential interaction is that, if initially there is a zero probability for the particle to be at a; > 0, as soon 
as < = 0"*", there is instantaneously, a finite, though very small, probability to find the particle at any point a; > 0. 
This non-local behavior of the Schrdingcr solution is due to its nonrelativistic nature and not a result of the quantum 
shutter setup 12 1. Application of the Klein— Gordon equation to the shutter problem [l| shows that the probability 



density is restricted to the accessible region x < ct {c is the speed of hght). 

Gerasimov and Kazarnovskii Q confined the initial wave in a finite region by introducing a second shutter at the 
point X = a. Godoy pointed out the analogy with FraunhofFer diffraction in the case of small box (compared to 
the de Broglie length), and Fresnel diffraction, for larger confinements. 

On the other hand, the phenomenon of quantum superarrivals for the reflection and transmission of a Gaussian 
wave packet in the case of a rectangular potential barrier while it is perturbed by either reducing or increasing its 
height has also been considered [8|-|l0| . There exists a finite time interval during which the probability of refiection is 
larger (superarrivals) when the barrier is lowered compared to the unperturbed case Q. Similarly, during a certain 
interval of time, the probability of transmission when the barrier is raised exceeds that for free propagation Q. Such 
studies suggest that the time-varying boundary conditions can give rise to interesting action-at-a-distance effects in 
quantum mechanics. 

Exactly solvable models give important points about the quantum many-body physics. The Thonks- Girardeau 
(TG) gas [ll] is a such model, which corresponds to infinitely repulsive bosons in one dimension (ID). This model is 
exactly solvable via a Fermi-Bose mapping, which relates the TG gas to a system of noninteracting spinless fermions 
in ID [13. 

In this context, by considering the problem of a confined particle in a linear potential, the aim of the present paper 
is to probe some aspects of the time-dependent boundary conditions that have remained hitherto unnoticed. 

The organization of the paper is as follows: In Section |ll] we review the solution of the Schrodinger equation for a 
confined particle in half space < x with the potential V{x) = Kx and consider the various time-dependent boundary 
condition of the problem. Section lllll brieflv deals with quantum statistics. Numerical calculations are presented in 
section HVl Section fVl contains concluding remarks. 
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II. A CONFINED PARTICLE IN HALF SPACE 0<x WITH THE POTENTIAL V{x) = Kx 

A. Stationay states 

The Schrodinger equation of the n*'' state of the particle wave function is 

(f(f>„(x) 27TI . X , / X „ 

—^J-^ + ^{E„-Kx)(l)n{x) - 0. (1) 

Performing the change of variable 



a 



f2mK\ 



1/3 



(2) 



we get the Airy equation, 



whit solutions, 



^+e0n = O, (3) 



cl)^{x) = N^Ai{-0 ■ (4) 
in which Nn is the normalization constant. The energy levels are determined by the condition (/),i(0) = 

En 



Ai 



-a- 
K 



= . (5) 



Then, the i?„ are determined by the zeros a„ of the Airy function 



En = -an— = -an — . (6) 

a \ 2m ) 

Thus, the solution of Eq. ([T]) is given by [isj . 

(/)„(a;) = V"-r77 — r^i(aa; + a„) , (7) 
Ai \a,n ) 

in which ^i'(a„) = ^^^^|x=a„- Now, we study the evolution of Eq. ([7]) under time-dependent Schrodinger equation 
in the following cases: 

B. Sudden removal of the wall and switching on the linear potential Vix) — Kx at a:: < 

Suddenly the wall at a; = is removed at time equal to zero, and a linear potential of the form V{x) ~ Kx is 
switched on in the region a: < 0. Now, the particle is in a linear potential in whole space. At any instant t > 
wavcfunction is obtained by, 

/•OO 

ip{x,t) = / dx' G{x,t\x',0)Mx') , (8) 



in which "00(2^) is the initial wavcfunction and is chosen as Eq. ([7]) and G(a;,t|x',0) is the propagator for the Hamil- 
tonian H = /2m + Kx and is given by p^ . 



G{x,t\x',0) = e^^(^-^')'-W(-+-')-^ . (9) 
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Putting Eq. dH) into Eq. one finds 



^nix,t) = Va-T-yj re ^ 2« 



yli'(a„) V '^TTiht Jq 

Inverse Fourier transform of tlie function Ai{k) ^ e«fc^/3 jg ^j^q Airy function, 



27r 



dk e"=^Ai(fc) 



(10) 



(11) 



So, 



^„(x,i) = = ^^^^ e-W-^ / d/fc e^fcV3 e''^-" / dx' e^^^-'^" 



271 Ai'{a„ 



2Triht 



2h ^ ^iakx 



2mht 







1 1 

A-Ky/a Ai'{ari 

in which X ~ x + Kt^ /2m — hkt/m and erfc(a;) is the complementary error function 



/ dk e"= e'''"'^/" e*S(^--^ 'erfc 



i + 1 [rn 



2 \ U 



X 



(12) 



C. Sudden removal of the wall and the potential (Free propagation) 

By sudden removal of the wall and turning off the potential at time equal to zero, the particle is allowed to propagate 
freely. Using Eq. ([9]) for JiT = 0, one finds 



1 1 

AiTy/a Ai'{an) j_ 



^ik^/3a='-iEt/n+ik(a„ /a+x) ^^^^ 



» + 1 [rn ' 



(13) 



in which, E = h^k^/2m and v = fik/m. From Eqs. ([T^ and ([T^ one gets the following relation between the density 
probability in the free and the interacting case: 



|^i^)(x,i)|^ = \^^^\x+'^,t)\^ 
Thanks to this relation we do numerical calculations only for the free case. 



(14) 



D. Sudden removal of the potential 



In this case the propagator is given by [17| . 

G{x,t\x',{)) ^ Gfrcc(x,i|a;',0)-Gft.eo(-x,i|a;',0) , 



(15) 



where Gircc{x^t\x' denotes the propagator of free particle and is given by Eq. © with K ^{). From Eq. (|13p . we 
have 



1 1 



A-K^ Ai'{i 
X ( e'''^erfc 



•n) J-oo 



i + 1 [rn 



2 V fit 



(x — vt) 



e '''^erfc 



i + 1 [rn 



2 V fit 



{~x — vt) 



(16) 



III. QUANTUM STATISTICS 



So far we have not consider quantum statistics. We shall next consider dynamics of the ground state of a system of 
N noninteracting spinless Fermions initially confined in half-space x > in the presence of a linear potential. Since 
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wavefunctions of noninteracting spinless fermions are antisymmetric under coordinate exchanges, their wavefunctions 
vanish automaticaUy whenever any Xj = Xk- 

tpFo{xi,...,XN,t) = -^det|^^:^2(i,i)V'«(a;j,<) • (17) 
The single particle density, normalized to N , is 

N 

p{x,t) ^ N j \i]Fo{,X,X2,.:,XN,t)\'^dX2—dXN='^\'lpn{x,t)\'^. (18) 



n=l 



It follows from the Fermi-Bose mapping theorem |18l-l2d| that the exact A^— bosons ground state wavefunction ipBo of 
a system of N bosons with hard-wall interaction is 

^pBo{xi,...,XN,t) = \lpFo{xi,...,XN,t)\ . (19) 



IV. NUMERICAL CALCULATIONS 



In this section we work in a unit system where fi — 1 and m = 0.5. We choose K = \ anywhere the potential is 
nonzero. Thus, from Eq. ^ the initial wavefunction is given by, 

V'o(a;) = (t>n{x) = —^ — -Ai{x + an) ■ (20) 



We use Kriiger's method [la, |21| to produce Airy functions. 



A. Sudden removal of the wall and the potential (Free propagation) 



Fig. [Ushows the probability density versus distance x for state n = 6 with energy eigenvalue Eq = 9.023 at different 
times. From this figure one can see that at long times probability density locates symmetrically around the origin 
and it spreads without motion. From Eqs. ^ and (j9]) for free propagation one has. 



and 




in which c.c. refers to the complex conjugate of the first term in the bracket. One can check numerically in the 
limit i — > oo. last line of the above expression is very small compared to the two first terms. As a result of which 
density probability would be an even function of x at long times. After a certain time t„, for even n, |'!/)„(a;, <)p has 
a local minimum at x = and its n maxima arc distributed symmetrically around a; = 0, but for odd n it has a 
local maximum at a; = and other n — 1 maxima locate symmetrically around the central maximum, i.e., (n — l)/2 
maxima at the right and (n — l)/2 maxima at the left of the origin. Another feature of [■(/'^(a;, t)P is that it has 
no node after a certain time of letting the particle to be free. From Eq. one can find this is also true for 

propagating of 0„ (x) in a linear potential. Another physically interesting quantity is the probability current density 
j = (h/ m)'^{ip* dip / dx) . It is important because arrival time distribution at a given location is given by its modulus 
sign in some approaches to this problem p^ . Fig. |3] shows the probability current density versus time at different 
locations. From this figure one can see j is negative for negative values of x and is positive for large positive values of 
X, but its sign changes with time for small values of positive x. Fig. [5] represents the evolution of the density profile 



5 



for a noninteracting spinless Fermioninc gas composed of iV = 6 atoms released from the trap. When the particle is 
trapped in the rt— th mode, and after being released, the probability density presents n maxima and by elapsing the 
time, they distribute symmetrically around the origin. As a result of which, for the case of A^— body system, single 
particle density presents N maxima, when particles are trapped and after being released and by passing the time 
these maxima locates symmetrically around the origin. 

B. Sudden removal of the potential 

Fig. [2]shows probability density versus distance x for state n = 6 at different times. It cen be seen that wavefunction 
spreads without motion. The height of the last maximum is larger than the n—1 first ones before turning the potential 
off, but by going time, the height of them goes to be the same except the last one, whose height is smaller than the 
others. \il)n{x,t)\'^ in Fig. [2] has always n nodes in spite of Fig. [1] This comes from the fact that in this case the 
wall is in place. Fig. |4] shows the probability current density versus time at different locations. From this figure one 
can see that j is positive for large values of x, but its sign changes with time for small values of x. Finally, Fig. [S] 
represents the evolution of the density profile for a Fermioninc gas composed of = 6 noninteracting spinless atoms 
released from the trap. Noting this figure, one sees the height of the maxima decreases by x for a given time. 

V. SUMARRY AND DISCUSSION 

In this paper we have considered numerical solution of the time-dependent Schrodinger equation for some time- 
dependent boundary conditions. Probability density, probability current density and dynamics of the ground state of 
a system of A^ noninteracting spinless Fcrmions were studied. When the particle is trapped in the n— th mode, and 
after being released, the probability density presents n maxima. As a result of which, for the case of A^— body system, 
single particle density presents A^ maxima, when particles are trapped and after being released. Since the stationary 
wavefunctions </)„ (x) given in Eq. ([7]) form a complete set of orthonormal eigenfunctions in the range (0 < x) , we can 
expand an initial wave packet ■!/'(x,0) which has an arbitrary shape for (0 < x) in terms of (j)n{x) 

oo 

(x) , (23) 

71=0 

where the expansion coefficients c„ is given by: c„ = 0* (x)V'(x, 0) dx. We have find time evolution of eigenfunctions 
4>n(x) which were denoted by ipn(x, t). Thus, time evolution of Eq. ([^^ is given by 

oo 

1p{x,t) = J2^ni'n{x,t) . (24) 

Time-dependent probability density \ip{x, t)\'^ is a complicated linear combination of [ipnix, i)P and the corresponding 
interferences amplitudes tp^{x^t)ipn'ix,t). 
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FIG. 2: Probability density versus distance x for state n = 6 at times (a) t = 0, (b) t = 10, (c) t = 100, (d) t = 500, (e) 
t = 1000 and (f) t = 2000. (Fee propagation in half-space x > 0) 




FIG. 4: Probability current density versus time t for state n = 6 at points (a) x 
(Free propagation in half-space a; > 0) 



= 5, (b) X = 



10, (c) a; = 15 and (d) x = 20. 
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FIG. 5: Density plots of the probability density p — X^^^j^ \iljn{x,t)f for a noninteracting spinless Fermioninc gas composed of 
N ^ 6 particles versus distance x at times (a) t = 0, (hYt = 10, (c) t = 100, (d) t = 500, (e) t = 1000 and (f) t = 2000. (Free 
propagation) 
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FIG. 6: Density plots of the probability density p — X^^^j^ \'<pn{x,t)f for a noninteracting spinless Fermioninc gas composed of 
N = 6 particles versus distance x at times (a) t = 0, (b) t = 10, (c) t = 100, (d) t = 500, (e) t = 1000 and (f) t = 2000. (Free 
propagation in half-space a: > 0) 



